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ABSTRACT 

Binary evolution predicts a population of helium core [M < 0.5 Mq) white dwarfs (WDs) that are 
slowly accreting hydrogen-rich material from low mass main sequence or brown dwarf donors with 
orbital periods less than four hours. Four binaries are presently known in the Milky Way that will 
reach such a mass-transferring state in a few Gyr. Despite these predictions and observations of 
progenitor binaries, there are still no secure cases of helium core WDs among the mass-transferring 
cataclysmic variables (CVs). This led us to calculate the fate of He WDs once accretion begins at 
a rate M < 10"^" Mq yr~^ set by angular momentum losses. We show here that the cold He core 
temperatures (Tc < 10^ K) and low M result in 10^^ Mq of accumulated H-rich material at the 
onset of the thermonuclear runaway. Shara and collaborators noted that these large accumulated 
masses may lead to exceptionally long classical nova (CN) events. For a typical donor star of 0.2 Mq, 
such binaries will only yield a few hundred CNe, making these events rare amongst all CNe. We 
calculate the reheating of the accreting WD, allowing a comparison to the measured WD effective 
temperatures in quiescent dwarf novae and raising the possibility that WD seismology may be the 
best way to confirm the presence of a He WD. We also find that a very long (> 1000 yr) stable burning 
phase occurs after the CN outburst, potentially explaining enigmatic short orbital period supersoft 
sources like RX J0537-7034 (Porb = 3.5 hr) and IE 0035.4-7230 (Porb = 4.1 hr). 

Subject headings: accretion, accretion disks — binaries: close — nuclear reactions, nucleosynthesis, 
abundances — novae, cataclysmic variables — stars: dwarf novae — white dwarfs 



1. INTRODUCTION 

For short enough orbital periods, the evolution of 
a low-mass star in a binary system up the red gi- 
ant branch can be halted by the onset of catastrophic 
filling of its Roche lobe. This exposes the low-mass 
(M < 0.5 Mq; e.g., 'Dominguez et al.lll999t iSalaris eFaLl 
I2002t iSerenein &: Fukugita ,2001) helium core of the 
red giant and trigge rs a common envelope event (see 
iSandquist et al.l[2000l for an overview) that brings the re- 
sulting cooling He white dwarf (WD) much closer to its 



stellar com panion (dc K oo]||1992tlde Kool fc R ittei"1993' 
llben fc Tut ukov 1993). The discoverv by iMarsh et al.i 
(|1995D that low-mass WDs have stellar companions sup- 
ports this origin for many He core WDs. 

Angular momentum losses (either from gravitational 
waves or stellar wind braking) after the common enve- 
lope event will drive the binary closer, with some reach- 
ing contact i n a few Gyr to become cat aclysmic vari- 
ables (CVs). iSchreiber fc Gansickd (|2003f) studied such 
pre-CVs, noting that a few have He WDs. The re- 
cently discovered system WD 0137-349 is a new pre-CV 
example with a 0.39 Mq WD and a 0.053 Mq brown 
dwarf com panion that will come into contact in 1.4 Gyr 
Maxted et al. 2006). Populat i on synthesis calcula tions 
de Kooll 11993 iPolitanol Il996l: IHowell et all [2Ci0ll) pre- 
dict that as many as « 20% of CVs with orbital pe- 
riods Porb < 2 hr should have He WDs. However, 
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where WD masses are inferred in CVs (|Pattersonll2001l : 
Townsley & Gansickc 2009), no such low masses have 
been identified. 

These considerations led iShara et al.l (|1993D to exam- 
ine classical nova (CN) events in these unusual systems. 
They calculated the outcomes for a 0.4 Mq He WD ac- 
creting at a rate of M = 10~^ Mq yr~^, finding accu- 
mulated masses of Mjgn = 9 x 10~^ Mq for a WD core 
temperature of Tc = 10^ K and a time between outbursts 
of 10^ yr. They evolved the WD for nine flashes, finding 
slow novae followed by a prolonged bright (L > 10'^ Lq) 
phase after the outburst, as the residual « 10""* Mq hy- 
drogen shell burns stably. The combination of a large 
Mjgn in a tight binary may m ake these events excep- 
tional. In the case of C/0 WDs. llben fc Tutukovl ([1993 ) 
suggested that such events would trigger a common en- 
velope event of some significance, potentially explaining 
the very bright "red variable " seen in M31 (|Rich et al.l 
[TgsllBoschi fc Munaril[200l . 

The known pre-CVs with He WDs provide excellent 
starting points for our work, so we begin in !j2] by not- 
ing the stability of mass transfer once contact is reached 
and deriving the subsequent M{t) historie s. We find 
that l ower accretion rates than studied by IShara et al.l 
(|1993f l need to be considered, so in ij3]we perform (and 
analyze) a few Gyr of time evolution at constant low 
M = 10-1° Mq yr-i and 10"" Mq yr"! on a 0.4 Mq He 
WD. Given the low A/, there is time for di ffusion to occur 
betw e en the fresh mate rial and the WD (jPaauette et al.l 
[l986t llben et al.l fl99l . which we discuss in the Ap- 
pendix. 

In f|4l we perform calculations of two specific scenar- 
ios motivated by the known He WD pre-CVs. Our focus 
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there is on the long-term evolution of their CN outbursts 
as well as the Tc evolution. We also discuss the WD 
surface temperatures when in quiescence as dwarf nova 
systems and the potential discove ry of He WDs through 
their pulsations (| Arras et al. !'2006). In fj5l we explore de- 
tails of the individual novae, describing the physics of the 
convective burning phase in H5.ll and the composition of 
the nova ejecta in H5.2[ which is not enhanced in CNO 
nuclei due to th e lack of a reservoir of underlying C/0 a s 
in normal CNe (jStarrfield et al.lll97alGehrz et al.lll998n . 
The CN event leaves a remnant burning envelope on the 
surface of the WD, allowing for a prolonged stable burn- 
ing phase (jShara et al.lll993[ ) that we detail in ^5.'S[ which 
may explain the few Porb ^ 4 hr long-lived supersoft 
sources (SSSs). We conclude in |j6]by highlighting where 
future CV and CN observations may reveal the long-lost 
helium core WDs in mass transferring binaries. 

2. STABLE MASS TRANSFER RATES 

The common envelope event that exposes the He WD 
leaves a wide range of WD and donor masses (Md) , only a 
fraction of which will become stable mass-t ransferring bi- 
naries once the donor comes into c ontact (|de Koo]||1992l : 
IPolitanol [1991 IHowell et all [200lf) . Table [J shows the 
parameters of those pre-CVs with probable He WDs 
whose companions will overflow their Roche lobes in 
iin < 20 Gyr, under the assumption of angular momen- 
tum losses due to gravitational waves. The prevalence 
of low-mass {Md < 0.25 M©) main sequence compan- 
ions is likely due to their slower inspiral under gravita- 
tional wave losses compared to those with more massive 
{Md > 0.25 M(7i) companions that would have stellar 
wind braking (jSchreiber fc Gansickd [20031 ) . As we will 
show in the following sections, all binaries in Table [T] 
will likely transfer matter stably after contact, leading 
to the accretion of cosmic-mix material onto a He WD 
at M = 1 - 4 X 10"" Mq yr^i. 

2.1. Stability of Mass Transfer 

When the main sequence star comes into contact, 
the possibility exists for dynamically unstable mass 
transfer, which occurs when the donor's Roche ra- 
dius shrinks faster or expands slower t han the donor's 
radial response to adiabat ic changes (jWebbinkl 119851 : 
IHiellming &: WebbinkI IT987I) . Low -mass mam sequence 
stars have small radiative cores and large convective en- 
velopes that make up > 50% of the star's mass, and 
stars with Md < 0.3 Mq are nearly fully convective. 
The adiabatic response of the radius of a fully convec- 
tive star to mass loss is d\nRd/d\nMd\g = —1/3, and 
using the Paczvhski (1967) approximation to the donor's 
Roche radius yields dlnR^/dluMd ~ —5/3 + 2q, where 
q = Md/M. Thus, dynamically unstable mass transfer 
occurs if the donor is fully convective and has a mass ra- 
tio Md/M > 2/3. For the maximum mass He WD with 
M ~ 0.5 Mq, this restricts the main sequence donor 
mass to Md < 0.33 Mq in order to avoid dynamical mass 
transfer. If the more exact Roche radius formula from 
lEggletorJ (|1983f ) is used, the limit changes only slightly 
to Md/M < 0.63. Binaries just at the limit of stability 
would have donors that fill their Roche lobe s as main se- 
quence stars {Md = 0.33 A/©, Rd — 0.33 i^©: lRibas et al.l 
[2OO81) with Porb = 3 hr. Any bloating of such a donor 



could drive Porb slightly longer, potentially to the 3.5 — 4 
hr seen in a few supersoft sources, as we discuss in 

All the pre-CVs in Table [1] have Md/M < 0.63 and 
will become stable mass transfer systems when they come 
into contact a time tin from now. We derive their mass 
transfer rates in i j2.2l There are certainly He WD post- 
common envelope systems with less extreme mass ratios 
that would trigger a dynamical mass transfer event at 
the onset of Roche lobe filling. These yield a red giant 
configuration, as the supply of freshly accreted H leads 
to a rejuvenation of the H burning shell typical of a He 
core on the first ascent of the red giant branch. 

2.2. Mass Transfer Rates for Stable Scenarios 

Our focus here is on those systems that undergo sta- 
ble mass transfer, with the WD accumulating mass over 
time and then ejecting most of it during the CN. We start 
with a total mass Mt = Md + M and orbital angular mo- 
mentum J — MdM{Ga/ MtY^"^ for an orbital separation 
a. Mass transfer is driven by the rate of orbital angular 
momentum loss, J, via 
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Since CNe lead to ejection of material from the binary, 
we must account for mass loss. We define an excavation 
factor / by setting Mt = fMd, so that f — 1 means the 
WD, on average, keeps a constant mass (i.e., the CNe 
eject the amount of matter that is accreted), and f = 
means that the WD keeps all the accreted mass. The 
resulting relation is then 
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Presumin g that the donor a lways fills its Roche lobe, and 
using the iPaczvhskil (|1967f ) formulation, we obtain 
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(3) 



We need the response of the donor star's radius to mass 
loss, Crf = d\nRd/d\nMd, which initially occurs on a 
timescale longer than the Kelvin-Helmholtz time, allow- 
ing the star to stay close to the main sequence (e.g., 
Cd ~ !)• However, the increase of the Kelvin-Helmholtz 
time as Md decreases leads to a nearly ad iabatic response 
(Cd ~ — 1 /3) late in the binary evolution (jKolb fc Baraffd 
[l999h . 

The resulting mass transfer rate is then 
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f_Md 
3 Mt 



(4) 



where we use the work of iKolb fc Baraffd (fT999l) to ob- 
tain Rd{Md) and Cd- Their Rd{Md) relation was cal- 
culated self-consistently for a, M = 0.6 Mq WD with 
gravitational wave losses alone. We also presume that 
gravitational-wave losses set J and consider a few spe- 
cific scenarios in Figure [TJ The two solid lines labeled A 
and B are the two scenarios we will discuss at length in 
521 Scenario A begins with a Md = 0.2 Mq donor, while 
scenario B begins with a Md = 0.05 Mq donor. Both 
assume that all matter accreted onto the WD is ejected 
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TABLE 1 

Helium WD Pre-CVs with tin < 20 Gyr (as of July 2009) 



System 


M 


Ma 


tin 


Ref. 




(Me) 


(Mq) 


(Gyr) 




SDSS J1435+3733 


0.48 - 0.53 


0.19 - 0.25 


0.75 - 1.1 


1,2 


HR Cam (GD 448) 


0.41 ±0.01 


0.096 ± 0.004 


1.5 - 2.1 


3 


WD 0137-349 


0.39 ±0.035 


0.053 ± 0.006 


1.7 - 2.5 


4 


SDSS J1529+0020 


0.40 ±0.04 


0.25 ± 0.12 


0.54 - 6.3 


5 


SDSS JOllO+1326 


0.47 ±0.2 


0.255 - 0.38 


12 - 19 


2 


SDSS J1724+5620 


0.42 ±0.01 


0.25 - 0.38 


13 - 20 


5 


RR Cae 


0.44 ±0.022 


0.182 ± 0.013 


18 - 23 


6 


SDSS J1212-0123 


0.33 - 0.48 


0.26 ± 0.03 


16 - 26 


7 



1. ISteinfadt et al.l ^OOB). 2. IPvrzas et al.l ||2009|'). 3. IMaxted et al] 119981) . 4. IMaxted et all (20QB) . 5. IRebassa-Mansergaa et aTl 112001) . 6. 
IMaxted et al.l H2007l) .~ INebot Gomez-Moran et al.l H2009I1 



during the CN event. The resulting accretion rates are 
< 4 X 10""'^^ Mq yr~^ once mass transfer begins. The 
trends that we discuss throughout the rest of the paper 
for 0.2 Mq donors will apply to larger donors as well 
because the accretion rates determined by gravitational 
wave radiation for low-mass main-sequence donors differ 
by a factor of a few at most. For example, the accretion 
rate for the most massive stably mass-transferring sys- 
tem {Md = 0.33 Mq, M = 0.5 Mq) is only 3 times larger 
at a given orbital period than our fiducial example with 
Md = 0.2 Mq and M = 0.4 Mq. 

We neglect any potential variations of M on the CN re- 
currence time due to prolonged periods of "hibernation" 
(|Shara et al.lll986f l. We further note that the Rd{Md) re- 
lation we use assum es a donor with sola r meta llicity. As 
iStehle et all {1997') and|Kolb_&Baraffe (1999) note, the 
secondary's composition does have an effect on its struc- 
ture and thus on the secular evolution of the CV. We have 
also neglected the small correction needed to re-adjust 
the Rd{Md) relation for our specific M history, which is 
lower by a factor of « 1.5 than that in iKolb &: Baraff^ 
(|1999f) due to our 0.4 Mq WD (they modeled a 0.6 Mq 
WD). This would tend to slightly shorten the minimum 
orbital period, as the mass-transfer timescale at a given 
donor mass is longer in our case. However, given that the 
current calculations do not match the now observed mini- 
mum orbital period in either case (iGansicke et al. 2009), 
we will neglect the corrections due to the secondary's 
composition and M history. 

3. REHEATING OF SLOWLY ACCRETING 
HELIUM WHITE DWARFS 

The previous section showed that the relevant ac- 
cretion rate s are much lower than that considered by 
iShara et al.l ()1993D . In addition, by the time the known 
pre-CVs come into contact and initiate stable mass trans- 
fer, the He WD core will hav e cooled to w 2-6 x 10^ K 
(lAlthaus fc Bcnvcnuto" 19971 . Extensions of the work by 
iTownslev fc Bildsten (2004) on CN ignition masses to 
lower M, M and yields Mig„ « 10"^ Mq. This im- 
plies CN recurrence times of 10^ — 10* yr during the 
many gigayears of evolution under accretion. The long 
evolution time allows for significant thermal coupling be- 
tween the freshly accreted H envelope and the He WD 
core during the accumulation of fresh matter. For this 
reason, we must track Tc and its evolution under the ac- 
tion of accretion for a prolonged period of hundreds of 
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Fig. 1. — The solid line in the left panel shows M as a function 
of ^orb assuming that all mass accreted onto the 0.4 Mq WD 
is ejected by CNe (/ = 1) for an initial = 0.2 Mq donor 
(scenario A; e.g., SDSS J1529+0020). The dashed line assumes 
that all accreted matter stays on the WD (/ = 0). The upper solid 
and dashed lines in the right-hand panel show the corresponding 
M{t) profiles. The lower solid curve in the right hand panel shows 
M for a donor with initial = 0.05 Mq (e.g., WD 0137-349 ), 
assuming that CNe eject all the accreted material (/ = 1; scenario 
B). 

CN events. 

3.1. Analytic estimate of thermal coupling between the 
envelope and core 

Due to the long accumulation and evolution time, 
much of the He WD core participates in the thermal 
balance of the CN cycle. We begin by estimating the 
timescale for heat transport between two mass shells at 
radial coordinates ro and r ([Henvev fc L'Ecuveilll969l ). 

3 ITfKcp^^/^ Y 
16asBC [Jra V T-^ / J 

where cp is the specific heat, asB is the Stefan- 
Boltzmann radiation constant, and k is the opacity. 
Electron conduction determines heat transfer in the de- 
generate He core, so k = k'T^/p^, where k' « 6 x 
lO"'^ g cm""' K-^ from fitting to iCassisi et all (j2007l )'s 
opacities at p = 2 x 10^ g cm~^ and T = 6 x 10^ K. 
When the core is isothermal with a constant cp, the ther- 
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Fig. 2. — Thermal diffusion time from the outside of the core to 
interior mass points for M = 0.3 Mq (dashed lines), 0.4 Mq (solid 
lines), and 0.5 Mq (dotted lines) He WDs. The isothermal core 
temperatures are marked. 

mal diffusion time becomes 



WasBcTc 
■ 1.6 X 10^ yr 

6 X 10^ K 



6 X 10-'^ g cm-4 K- 



cp 



f r-rp 
iks/'imp) V109 cm 



, (6) 



wliere we fiave normaiized cp to the dense iiquid vaiue 
for pure helium (i.e., an internal energy of iksT per ion, 
where ks is Boltzmann's constant and is the baryon 
mass).^ 

We numerically integrate equation ([5]) with the full de- 
pendences of K and Cp for varying masses and isother- 
mal temperature to find the thermal diffusion times from 
the core-envelope boundary to a given mass coordinate 
within the WD. These are shown in Figure [5] for 0.3 Mq 
(dashed lines), 0.4 M© (solid lines), and 0.5 Mq (dotted 
lines) He WDs. Each set of lines shows different core 
temperatures. From top to bottom, these are 3 x 10^, 
10^, and 3 x 10^ K. The only dependence on M in equa- 
tion ([6]) is from the radius of the core-envelope interface, 
which are all within 10 — 20% of each other for these low- 
mass WDs. Our simple estimate in equation ^ is fairly 
close to the integrated results in Figure [21 a core tem- 
perature of 6 X 10^ K yields « 2 x lO'* yr, and raising 
Tc by a factor of 10 decreases the thermal diffusion time 
by almost a factor of 10, as predicted by equation ([6]). 

3.2. Numerical Evolution of the Thermal State for 
Constant M 

In this section, we describe several time-dependent cal- 
culations of constant M onto a 0.4 Mq He WD in order to 

^ For low-mass He WDs, the Coulomb coupling parameter F is 
10, so the ions are in a liquid state, but not necessarily deep into 

the 3kgT limit. Thus, the actual values for cp are between the 

ideal gas and liquid values. 



discuss the thermal evolutio n and compare to the work of 
iTownslev fc BildstenI ()2004D . The specific outcomes for 
the M histories appropriate to CVs calculated in ^2.2\ 
will be discussed in [JH 

Numerical calculations were performed with the 
hydrodynamic La gr angia n stellar evolution code of 
iPrialnik fc Kovet3 ()1995f ). which contains an ex- 
tended nuclear reaction network of 40 elements up 
to 3ip (ICaughlan fc Fowled [Tosl lAngulo et all [ToM 



iTimmesI I1999D an d implemen ts the OPAL opacities 
(jlglesias fc Rogersl 119931 119961 ). Convection is treated 
by mixing l ength theory with a mixing-length parame- 
ter a = 2 (IMihalad [19781 : IMazzitellil \l97§ . IVandenbCTgl 
[l98l . Inter -species diffusion is allowed to occur (see 
the Appendix for discussions of diffusion), but no 
other mixing processes are presumed. N o ^He is in- 
cluded in the newly acc reted material (iSharal 119801 : 
ITownslev fc Bildsteni[200llShen fc Bildstenll2009a[ ). but 
its presence in the incoming material should be negligi- 
ble as the donor is unevolved. The CN mass-loss phase 
which follows the phase of rapid expansi on uses a stead y 
optically-thick supersonic wind solution (|Prialniklll986f) . 
Common envelope effects are neglected. 

The 0.4 Mq He WD has a core composition of 98% He 
and 2% ^*N by mass and is surrounded with an initially 
thin envelope of solar composition. Three cases were 
studied, with different M and initial T^: M = 10"^° 
and 10~" Mq yr^^ for initial = 6 x 10*^ K, and M = 
10"" Mq yr-i for initial = 8 x 10^ K. We numerically 
follow the evolution through hundreds of nova outbursts, 
with ea ch nova cycle calcu lated through all evolutionary 
phases. lYaron et al.l (j2005f ) found that novae at these low 
accretion rates tend to erode the WD (/ > 1). Therefore, 
after each outburst and resulting envelope ejection, the 
outer mass shell of the stripped WD core is split into 
multiple shells of < 10~^ Mq so that the next nova cycle 
can be accurately calculated. 

Equation ([6]) and Figure [21 indicate that the timescale 
for thermal coupling between the core and envelope is 
of order the timescale between nova outbursts for low 
M = 10"" Mq yr~\ although not for the higher 
M = 10"^" Mq yr~^ case, which has a shorter accu- 
mulation time of « 6 X 10^ yr. However, because the 
envelope has a much lower mass (and thus a smaller ther- 
mal content) than the core, it will take many nova cycles 
to significantly increase Tc- The thermal content in the 
nearly ideal gas envelope of mass Mcnv with mean atomic 
weight /i w 0.6 is ii^th.onv = 3McnvkBT/2fj,mp « 4 x 
lO'^^ erg (T/IO^ K)(Mo„v/10~^ Mq), whereas the ther- 
mal content in the He core is -Eth.core ~ "HMkBTd 'imp k, 
5 x 10"*^ erg (TjlO^ K)(M/0.4 Mq). The 100-fold ratio 
of these thermal energies implies that it will take ~ 100 
nova cycles to significantly heat the core if the only en- 
ergy source is thermal. Some nuclear "simmering" could 
shorten the required time, as we discuss later in this sec- 
tion. 

Figure [H shows the time evolution of Tc for the 
three cases, all eventually reaching a final equilibrium 
core temperature. The core temperatures have been 
smoothed in order to make the evolution more clear. The 
dashed line is for M = 10"^° Mq yr~^, which reaches 
an equilibrimn temperature of Tc ~ 9 x 10^ K and sur- 
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3x10' 



4x10'' 



1 2 3 4 5 

Time Since Accretion Started (Gyr) 

Fig. 3. — Time evolution of the core temperature, T^, of the 
0.4 Mq He WD for two different constant accretion rates and 
initial core temperatures. The solid (dotted) lines shows the 
M = 10~^^ Mq yr~^ case with initial core temperatures of 
8 X 10^ K (6 X 10® K). The dashed line is for an initial tempera- 
ture of 6 X 10® K, but M = 10~l° Mq yr~^. Lines are smoothed 
on a timescale of several nova cycles in order to make the secular 
evolution more clear. 
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Fig. 4. — Time evolution of the quiescent surface luminosity, Lq, 
for a 0.4 Mq He WD accreting at 10~^^ Mq yr"^ after it has 
reached its equilibrium core temperature. 



face luminosity in disk quiescence Lq k, 7 x 10^'^ Lq. 
The time needed to reach this equiUbrium is about that 
expected from the earher considerations regarding the 
thermal content and heating of the core from the accu- 
mulated envelope, as « 0.2 Mq of matter was put on the 
WD through this evolution, triggering 200 CNe. 

The evolutions for lower M = 10~^^ Mq yr~^ are 
shown by the solid and dotted lines in Figure [3] for two 
different initial core temperatures. For these cases, the 
initial Tc is higher than the equilibrium value; accretion 
serves to slow the cooling of the WD as opposed to heat- 
ing the core in the higher M = 10~^° Mq yr~^ case. The 
equilibrium core temperature is Tc = 5.5 x 10^ K with 
a quiescent surface luminosity of « 7 x 10"'' Lq for 
the majority of the accretion phase of each nova cycle, 
as shown in Figure HI which plots the evolution of the 
quiescent surface luminosity during a typical nova cycle 
after the equilibrium core temperature has been reached. 



iTownslev fc BildstenI (|2004f ) and lEpelstain et~all 
(|2007f ) showed that equilibrium core temperatures 
would be reached under the action of accretion and 
CN events but did not co nsider WD masses this low. 
ITownslev fc BildstenI (j2004[ ) showed that the equilibrium 
core temperatures could be understood as that temper- 
ature where the cooling of the core (mostly during the 
accumulation stage) would be matched by the heating of 
the core (typically late in the accumulation, approaching 
the CN event). They highlighted two energy sources 
during accumulation: "compressional" heating (really 
entropy loss as material is advecting into the star) and 
slow nuclear burnin g. Following the derivation in Ap- 
pendices A and B of ITownslev fc Bilds tcu (20 Q3), which 
in turn follows iN omoto fc Sugimotol (| 1977i) . iNomotol 
^982), and Her nan z et al.| (| 19881) . vields the equivalent 
of ITownslev fc Bildstcn (2004)'s equation (B4) for the 
"compressional" heating luminosity, 

L,«4.5^M 



',6.1 X 10"* Lr. 



M 



Iff K, 



10- 



Mq yr- 



,(7) 



after accounting for the differences in the Coulomb cou- 
pling parameter, mean atomic weight, and tempera- 
ture in the core, and the larger Mjgn. This predicts 

3 T 

Lq 

10 fl/f^ 



Lcam-n = 3.4 X li) ^ Lq and Lcomp = 5.5 X 10 



^0 

'1 



comp 

for M = 10"" Mq yr-1 and M = 10""' Mq yr-\ re- 
spectively, after their equilibrium core temperatures have 
been reached. These estimates for the outgoing lumi- 
nosity are « 50% lower than Figure 4 shows for M = 
10~^' Mq yr"^ and « 20% lower than the numerically 
calculated Lq ^ 7 x 10"^ Lq for M = lO"'" Mq yr"'. 

These differences in luminosity can be attributed to 
our numerical c alculation's inclusion of nuclear simmer- 
ing. lTownslev~fc Bildsten (2004) showed the potential of 
nuclear energy release during accumulation especially for 
massive envelopes like these (see their Figure 3). How- 
ever they did not carry out this specific calculation. Our 
time dependent calculations allow us to see how the nu- 
clear energy source can come into play. 

Figure [H] shows snapshots of the thermal profile dur- 
ing different stages of the 10* yr CN cycle for M = 
10"^' Mq yr"^ after the WD has reached its equilibrium 
core temperature. The three dashed lines (top to bot- 
tom) are for times 5 x 10^ yr, 7 x 10^ yr and 10^ yr, after 
the CN. The temperature inversion shows that heat flows 
from the hot post-CN envelope into the region beneath 
the burning layer during this time, to exit later during 
the first part of the accumulation phase. After this heat 
has been radiated away, the thermal profiles (three lower 
solid lines at 3 x 10^ yr, 5 x 10^ yr and 7 x 10^ yr after the 
outburst) begin to resemble cooling WD profiles. Accu- 
mulation of material (and nuclear simmering) redevelops 
the temperature inversion (top solid line at 9 x 10^ yr), 
and the dotted line shows the profile at 9.95 x 10^ yr, 
when convection has initiated. Simply put, the excur- 
sions of the temperature during the CN event and the 
long accumulation time allow for some nuclear energy to 
be transported into the outer layers of the WD core over 
time, to be radiated later in the accumulation stage. 
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Fig. 5. — Density and temperature evolution during the period 
after a CN and the subsequent accumulation to the next explosion 
for a 0.4 Mq WD accreting at M = IQ-li Mq yr'^. The three 
dashed lines are for times 5 X 10^ yr, 7 X 10® yr, and lO'^ yr after 
the CN, from top to bottom. The solid lines, from bottom to top, 
are for times 3 X lO'' yr, 5 X 10^ yr, 7 X lO'' yr, and 9 X 10^ yr 
after the last outburst. The dotted line is the profile at a time of 
9.95 X 10^ yr, when convection has started. Bullets show the base 
of the accreted envelope. 

4. THERMAL EVOLUTION FOR CV SCENARIOS 
AND DWARF NOVAE IN QUIESCENCE 

As discussed in [21 stable mass transfer for low mass He 
WDs requires a low-mass companion that will transfer 
mass at the rate set by gravitational radiation losses. 
The more rapid mass transfer systems (i.e., above the 
period gap) that have stellar wind braking will be rarer. 
Our focus in this section is on the evolution of the two 
scenarios discussed in building on the intuition gained 
in gi 

4.1. Thermal evolution 

Figure [H shows the mass accreted between nova out- 
bursts (solid lines), envelope mass above the position 
where convection first begins during the thermonuclear 
runaway (dotted lines), and the ejected mass (dashed 
lines) for scenarios A and B of Figure [T] Even though 
the mass transfer rates and core temperatures differ in 
these scenarios, the typical CN ignition masses change 
very little. Most importantly, diffusion of the H into the 
He and the resulting convection during the outburst lead 
to an ejected mass that is slightly larger than what was 
accreted. Hence, our earlier approximation that the ac- 
creted mass is ejected was appropriate. 

The thermal evolution of the WD core is more complex 
due to the changing M in these realistic scenarios, as 
shown in Figure [7l The top panel shows the Tc evolution 
for which the initial donor is a 0.2 Mq main sequence star 
(scenario A; solid line), which exhibits an initial heating 
to a near-equilibrium core temperature of Tc = 6.8 x 10^ 
K and a luminosity of w 2 x 10~^ Lq (solid line in 
bottom panel) during the 4 Gyr when the accretion rate 
is nearly constant at M = 3 x 10~^^ Mq yr~^; this value 
of Lq is in agreement with the prediction of Equation ^ . 
This case is intermediate to the two calculations with 
constant M discussed in During these 4 Gyr, the 
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Fig. 6. — Accreted mass (solid lines), envelope mass above the 
position where convection first begins (dotted lines), and ejected 
mass (dashed lines) for the two binary evolution cases of Figure 
[T] Due to the diffusion of H into the He during accumulation, 
the initial convective envelope mass and ejected mass are always 
slightly larger than the amount of mass accreted since the last 
outburst. 
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Fig. 7. — The top panel shows the mass-averaged WD core tem- 
perature as a function of time for a 0.4 Mq He WD with initial 
Tc = 6 X 10® K for the two M histories (scenarios A and B) of Figure 
[1] The solid (dashed) line is for initial = 0.2 (0.05) Mq donor 
stars. The middle panel shows the time between CN outbursts for 
the same evolutionary scenarios. The bottom panel shows the evo- 
lution of the WD's surface luminosity in disk quiescence, Lq. For 
this plot, we use the minimum of Lq during each nova cycle, which 
because of the broad minimum is indicative of the 'typical' value 
(see Figure |4]l. 



CN recurrence time (solid line in middle panel) is 3 x 10^ 
years for about 130 CN events. After 4 Gyr, the donor 
begins to expand under further mass loss, and M declines 
as the binary passes through the period minimum. The 
WD then cools, and the recurrence time becomes quite 
long. 

The constantly declining M of scenario B (where the 
initial donor is a 0.05 Mq brown dwarf) gives a different 
thermal outcome, as the WD constantly cools (dashed 
line in upper panel) and its luminosity declines (dashed 
line in bottom panel). There are only 25 CNe in 6 Gyr; 
indeed, the donor itself only has enough mass for about 
50 CN events. 
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4.2. Dwarf novae in quiescence 

Both of these scenarios have low enough M that the 
accretion disk should be unstable, resulting in dwarf nova 
outbursts ()ShafteHll992[ ) and long periods of quiescence. 
When the accretion disk is in outburst, the accretion 
luminosity 

GMM 



2R 
i 0.004^0 



M 



0.4 Mq ^ 
M 



10 cm 
R 



10 



-11 



Mq yr-1 



(8) 



is an order of magnitude larger than the luminosity in 
the WD envelope, because the gravitational specific en- 
ergy, GM/R, is much larger than the thermal specific 
energy, kT^/ firrip, which sets the scale of the envelope's 
thermal profile. However, when the disk is in quiescence, 
the excess Lace is quickly radiated away on the thermal 
timescale near the WD's photosphere, and the luminosity 
of the system is dominated by the WD's quiescent sur- 
face luminosity, Lg, which origina tes from deep beneath 
the photosph ere (iPiro et al.ll2005D . 

Ever since iSionI ( 1999?) discussed the value of Toff 
measurements of dwarf novae in quiescence, observers 
and theorists have vigorously pursued this impor- 
tant diagnostic. Most recently, .Townslev fc Gansickei 
(|2009r ) summarized the o bservations and improved 
the ear lier theoretical work (iTow nslev fc Bildsten' ^20021: 
Godon & Sion 200l ITownslev fc Bildstcn 2003, 2005; 



Piro et al.ll2005i) . In FigureO we have added the Toff evo- 



lutions for our two sc e nario s to the original figure from 
ITownslev fc Gansickd (|2009[ ). Some data are marginally 
consistent with the solid line (scenario A), whereas no 
observed WDs are as cold as predicted from our scenario 
B; however, it is unclear if the absence of these cold sys- 
tems is physical or due to selection effects. It would be 
of interest to learn whether those cold systems near the 
solid line have any other evidence for having low-mass 
accretors. 

The continual discoveries of CV s in quiesceiice by 
the Sloan Digi t al Sky Survey fe.g.. iSzkodv et al.l 120091 : 
IGansicke etall '2009*) and, in the near future, by 
SkyMapper ([Murphy et al. 2008) will certainly provide 
new op portunities to reveal a ccreting He WDs. In par- 
ticular, IGansicke et al.l ()2009f )'s recent discovery of the 
expected "pile-up" of CVs at the 80 — 86 min orbital 
period minimum has alleviated long-standing concerns 
regarding binary evolution. Hence , as we noted in the i n- 
troduction, « 20% of the CVs in Gansick e etal\ (|2009f l's 
compilation should harbor a He WD. T he discovery and 
study of accreting WD pulsators (see iMukadam et al.l 
|20M for an updated list) may well be our best hope, 
as lArras et al.l (|2006D noted that low-mass WDs have 
colder effective temperatures for pulsation than higher 
mass WDs, and our Toff calculations shown in Figure 
5| are slight l y cold er than the blue edge calculated by 
' Arras et al.l (|20061 for low gravity (e.g., low-mass WDs). 
Thus, a prevalence of pulsators at low Tgs may be an 
indicator of He WDs. Certainly more work is needed to 
make this connection clear, but the rapid increase in the 
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Fig. 8. — Adopted figure from ITownslev fc Gansicki 1)20091 ). 
Solid (triangles) points are Tgff measurements from non- 
magnetic (magnetic) WDs in CVs below the period gap from 
ri'ownslev fc Gansickei 120091) . S haded regio ns are the predi cted 
effective temperatures found by ITownslev fc Gansick^ 1200911 for 
0.6 Mq and 1.0 Mq WDs. The solid and dashed lines are our 
calculations for a 0.4 Mq He WD with a 0.2 Mq and 0.05 Mq 
donor, respectively. We use the minimum in the quiescent surface 
luminosity for calculating T^f{. Figure |4] shows that for 90% of the 
CN cylo, the value of T^ff is no more than 15% higher than the 
value we plot. 

discovery of such pulsators is bound to reveal a few new 
systems worthy of intensive study. 

5. PHYSICS OF CLASSICAL NOVA OUTBURSTS 

In the previous sections, we studied the secular evolu- 
tion of the WD. We now turn our attention to details of 
the individual nova outbursts, focusing on the evolution 
of the convective burning phase, the composition of the 
nova ejecta, and the WD's post-nova appearance as a 
supersoft x-ray source. 



5.1. Evolution of the convective envel 
maximum temperatures 

For most of the accretion phase of the nova cycle, the 
thermal profile is set by "compressional" heating (see 
H3.'2\i . However, once the base of the envelope becomes 
dense and hot enough, the energy generation rate from 
nuclear burning becomes large enough that radiative dif- 
fusion can no longer effectively transport the luminosity. 
The radiative envelope is then transformed into a convec- 
tive zone whose entropy is increased by continued nuclear 
burning until mass is lost via a radiatively-driven wind, 
and the nova outburst is observed. We now give a brief 
overview of the relevant physics during this stage of the 
CN evolution; for a more complete explanation of simi- 
lar physics in He -burning convective envelopes, see §3 of 
IShen fc Bildstenl (|2009b) . 

Figure ini shows an example of the evolution of a convec- 
tive envelope as it is heated by nuclear burning, numeri- 
cally calculated for a 0.4 Mq WD that accreted mass at 
a rate of lO"^'^ Mq yr~^ after it has reached its equilib- 
rium core temperature of 5.5 x 10^ K. Dotted lines show 
temperature-density profiles of the envelope at 4.5 x 10®, 
5900, and yr before the convective zone reaches the top 
of our numerical grid, and solid lines show profiles 0.04 
and 2.1 yr after this time, with more evolved envelopes 
at lower densities. The dashed line is the adiabatic slope 
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Fig. 9. — Density and temperature evolution during a CN on 
a 0.4 Mq He WD. The dots display the evolution in time (from 
right to left) of the density and temperature at the base of the 
convective zone during the outburst. The uppermost dotted line 
is the profile when the convective zone reaches the surface of our 
numerical grid. The dotted lines beneath that are for times 5900 yr 
and 4.5 X 10® yr prior. The solid lines show the profile at times 0.04 
(when the base temperature reaches its maximum) and 2.1 yr after 
convection reaches the surface. The dashed line is the adiabatic 
slope, 9 In r/9 In p| J,, which is, as expected, nearly parallel to the 
convective region. The radial expansion of the outermost shell 
reduces the pressure at the base of the envelope, allowing for the 
underlying helium to adiabatically expand and cool (as is clear 
from the two solid lines at densities above 10'^ g cm~^. 

for an ideal gas, showing that the convective envefopes 
are indeed isentropic. 

During the initial heating of the convective zone, the 
envelope is geometrically thin, with a scale height h — 
Pb/pbd niuch less than the WD radius i?, and so the 
pressure at the base remains nearly constant: Pb « 
G'MMenv/47ri?^. Thus, as entropy is injected into the 
envelope during this phase, the temperature at the base 
increases and the density decreases while maintaining a 
constant pressure. Bullets in Figure [5] demarcate the 
base of the convective zone, which, as predicted, initially 
becomes less dense and hotter at constant pressure as 
the entropy increases. 

Further entropy injection causes the scale height to 
grow until it becomes a significant fraction of the WD's 
radius, and the envelope transitions to being geometri- 
cally thick. This leads to radial expansion of the en- 
velope and a subsequent decrease in the base temper- 
ature. Thus, there is a maximum base temperature 
for the convective envelope, which dep ends only on M, 
Menvi and the envelope's composition (jShen fc BildsterJ 
12009b'). For the model shown in Figure [SI the maximum 
temperature is « 10* K. 

Figure [TUl shows analytic calculations of the maximum 
base temperature as a function of the convective enve- 
lope's mass for a 0.4 Mq He core. The four lines are 
calculated with different hydrogen mass fractions, X, in 
order to parametrize the amount of core-envelope mix- 
ing. Increasing the helium abundance and thus the mean 
atomic weight /i decreases the number density of parti- 
cles, and so T is higher for a given pressure. Doubling 
X decreases the maximum base temperature by w 15%. 
Note that the temperatures for these models never reach 
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Fig. 10. — Maximum base temperatures vs. envelope mass for 
convective envelopes on a 0.4 Mq He WD. The curves are for 
different values of the hydrogen mass fraction, X, in the convective 
envelope, which depend on the amount of core-envelope mixing. 



TABLE 2 

CN EjECTA Mass Fractions on He and C/0 WDs 



Element 


Solar 


Ejecta (He WD) 


Ejecta (C/0 WD) 


H 


0.70 


0.62 


0.59 


He 


0.28 


0.36 


0.30 


Z 


0.017 


0.021 


0.11 


12c 


3.9 X 10"^ 


0.27 X 10"^ 


2.3 X 10"^ 


13 c 


0.43 X 10""* 


0.87 X 10-* 


7.7 X 10-"' 




0.10 X 10^2 


0.87 X 10-2 


5.7 X 10-2 




3.6 X lO"'* 


0.28 X 10"® 


2.4 X 10"® 


16 o 


0.94 X 10^2 


1.0 X 10^2 


4.8 X 10^2 




3.5 X 10-® 


1.4 X 10"=* 


3.0 X 10"^ 



values high enough to ignite the underlying He WD. Such 
an outcome would be exciting as it would likely lead to 
the birth of a He-burning star in a tight orbit. 

5.2. Composition of the nova ejecta 

Since the maximum base temperature in the convective 
envelope is not very high (w 10* K), the main channels 
for nuclear burning during the outburst are the p-p chain 
and CNO cycle. While CNe on C/0 WDs can dredge up 
large amounts of ^^C and ^^O nuclei from the underlying 
WD core, the only source of CNO nuclei for CNe on 
He WDs is from the accreted material and any trace 
CNO left in the He core. As a result, the ejecta has a 
metallicity, Z, close to the solar value, with an excess of 
H and He and a deficiency of ^^C and ^^O as compared 
to CNe on C/0 WDs. 

Table [2] shows the mass fractions for solar composition 
and for the numerically calculated ejecta of CNe on a 
0.4 Mq He WD accreting at M = 10"" M© yr"! and a 
0.65 Mq c/0 WD accreting at M = 10"^ Mq yr-\ 
computed after many nova cycles. Both calculations 
yield ejecta composed of ~ 10% core material and ~ 90% 
accreted material. This mixing results from the initiation 
of convection below the original core-envelope interface 
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due to diffusion during the accretion phase. However, be- 
cause the cores are composed of different material, Z is 
essentially solar for CN ejecta on He cores with a CNO 
breakdown that roughly reflects the abundances of the 
equilibrium CNO cycle, whereas the ejecta from the C/0 
WD CN has a metallicity that is 5 times larger than solar 
with overabundances of all of the CNO elements. Note 
that the specific isotopic ratios depend strongly on the 
reaction rates that are used and should thus be regarded 
as a rough guide. However, the abundance trends of H, 
He, and metallicity are a robust result. 

Classical novae whose ejecta exhibit near-solar metal- 
licities, and enhanced He and CNO-processed abun- 
dances can thus help to identify these systems. Of the 
currently small number of ~ 25 CNe wi th measured 
ejecta abundances (for com pilations, see IGehrz et al.l 
Il998t 'Hachisu h Katol 120061 ). none definitively show 
these observational signatures. However, given the fact 
that there should be far fewer CNe from He WDs as com- 
pared to C/0 WDs (see ^for an estimate of the rates), 
this is not surprising. 

5.3. Post-nova supersoft phase 

During a nova outburst, mass is ejected from the con- 
vectively mixed envelope (containing both accreted and 
core material) via a radiatively-driven wind. This mass 
loss continues until the envelope re-establishes a hydro- 
static solution, with a luminosit y close to the red giant 
core mass- luminosity relation of IPaczvhskil (|1970f ). The 
H-rich envelope undergoes quasi-steady nuclear burning 
with a nearly constant luminosity, decreasing in mass as 
H is converted to He until the remaining envelope be- 
comes too small to maintain a steady-state burning so - 
lution (Tuchm an fc Truran|[T99llSala fc Hernanzll2005l ). 
For near-Eddington luminosity and a photospheric ra- 
dius roughly equal to the radius of a 0.4 Mq core, the 
effective ter nperature yields a bright UV or s upersoft x- 
ray source (|Kahabka fc van den Heuveilll997D that lasts 
for as long as the envelope burns stably. 

Specifying the core mass, envelope composition, and 
luminosity, and demanding hydrostatic equilibrium and 
continuity of the radius and pressure at the core-envelope 
interface uniquely determine Menv in t his quasi-steady 
burning phase (see , e.g., Figs. 8 a nd 9 of lSala fc Hernanj 
120051 and Fig. 2 of iNomoto et al.ll2007D . As the envelope 
evolves, the luminosity remains nearly constant, decreas- 
ing slowly while Tcff increases and the envelope mass de- 
creases until a minimum envelope mass is reached, which 
coincides with the maximum Tcff- 

This initial phase of evolution is shown as solid lines 
in the H-R diagram in Figure [Tl] for He WD masses 
of 0.2, 0.25, 0.3, 0.4, and 0.5 Mq. The evolution for 
these solutions begins to the right of the figure at the 
point where the photosphere is just inside the He WD's 
Roche radius assuming Md/M = 1/2 and Porb — 4 
hr, which we take to be an estimate of the end of the 
mass-loss phase. The WD evolves towards the maximum 
Teff an d minimum M e nv at t he left. The equa t ion of 
state (ISaumon et al.l 119951 : iTimme s & Swestyl (2003, 
Rogers (^r^ avfonovl 2002[). opa city (Iglcsias fc Rogers 
19931 11996 : .Ferguson et al.l [2005.: .Cassisi et al.i 
2007 ) , electron sc r eening (iGraboske et al.l 11971: 




5.0 4.8 4.6 
log Teff [K] 

Fig . 11. — Evolution in the H-R diagram during the quasi-steady 
burning phase following a nova outburst. Each solid line shows the 
relation between L and T^g for a stably burning model for He core 
masses of 0.2, 0.25, 0.3, 0.4, and 0.5 Mq. Solid lines are truncated 
on the low T(,jf side where the radius of the photosphere equals 
the Roche radius assuming a mass ratio of 1/2 and Porb = 4 hr. 
The dashed extensions are thermally unstable branches that are 
not physically realized. 



Alastuev fc Jancovicil 11978 : 'Itoh et al. J 119791) , nu- 



clear burning network ((Caughlan fc Fowled [19881: 



lAngulo et a l."1999': ^ Timmesl '1999'). and neutrino cooling 
(jltoh etliLl.1996.) for these models are calculated with 
the MESA code package® (Paxton et al., in prep.). The 
mass composition of the envelopes is 90% solar compo- 
sition (as defined by Lodders 2003), which is accreted 
from the donor star, and 10% He core material due to 
convective mixing during the thermonuclear runaway, 
chosen to roughly match the ejecta abundances in Table 
[21 The CNO elements are assumed to be in chemical 
equilibrium as determined by their CNO-cycle burning 
rates, as their reaction timescales are shorter than the 
envelope's evolutionary timescale during this post-nova 
phase. Table [3| shows various properties of the envelope 
at the maximum Teff and minimum Menv for several 
core masses, where rsss is defined below, as well as 
the envelope mass at Roche lobe filling, Menv.RL- The 
maximum Toff ranges from several to tens of eV. 

While solutions exist for lower luminosities than that 
at the maximum Teff and minimum Menv, these solutions 
are not physically relevant because they require the enve- 
lope m ass to increase as time progresses (jSala fc Hernand 
120051 ). Furthermore, a stability anal ysis shows that 
these solu tions arc thermally tmstablc ('Paczvnski Hl983l : 
INomoto et a l. 2007; Shcn & Bildsten 2007). These unre- 
alized solutions are shown as dashed lines in Figure 1111 
In reality, once the minimum envelope mass has been 
reached, nuclear burning declines, the envelope mass re- 
mains constant, and the WD evolves to the WD cooling 
track until accretion leads to the next nova outburst. 
Figure [T^ shows the time evolution of the integrated nu- 
clear luminosity in the envelope and the quiescent surface 
luminosity from our numeric calculation after a typical 
nova outburst. The quasi-constant luminosity phase lasts 
for ~ 1000 yr before burning declines. 

^ http://mesa.sourceforge.net/ 
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TABLE 3 

Post-Nova Conditions at Maximum T^s / Minimum Mt 



M 






L, 




TSSS 


{Me) 


(10-3 Mq) 


(10-3 Mq) 


(Le) 


(IQS K) 


(yr) 


0.2 


3.1 


2.8 


5.1 


0.26 


3.7 X 10^ 


0.25 


1.1 


0.92 


46 


0.51 


1.3 X lO'* 


0.3 


0.52 


0.37 


2.4 X 10^ 


0.85 


1.0 X 10^ 


0.4 


0.20 


0.10 


1.4 X 103 


1.6 


4.8 X 103 


0.5 


0.11 


0.044 


3.3 X 103 


2.3 


8.9 X 10^ 




1000 2000 
Time (yr) 
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Fig. 12. — Integrated nuclear luminosity (dashed line) and quies- 
cent surface luminosity, Lq, (solid line) versus time after a typical 
nova outburst on a 0.4 Mq He WD accreting at 10—^^ Mq 
after it has reached its equilibrium Tc. The constant luminosity 
plateau lasts for 1000 yr, after which the nuclear luminosity 
drops sharply and the WD cools. 



We define a characteristic time the evolved post-nova 
spends near its maximum Tcff as rsss = McnvEnuc/ 
where E^uc ~ 4.1 x 10^^ erg g"^ is the energy per 
mass obtained by converting H to He in material with 
X — 0.68. As shown in Table [3l these timescales are 
as long as 5000 yr for a 0.4 M© core. Given this long 
duration, it is of interest to look for long-lasting (> 10 
yr) supersoft sources with short orbital periods < 4 hr. 
Possible systems like this are RX J0537.7-7034, which 
has Porb = 3.5 h r and a dynamic mass estimate of 
M = 0.4-0.8 Mq (lOrio fc QgelmanlfTOOa iGreiner et alj 



Orio et alj 



2000i), an d IE 0035.4-7230 ( P„rh = 4.1 hr; . 
1994; Sch midtke eFall 119961 ; Ivan Teeselinget al. 1998), 
whose closest non-LTE spe ctral fit y ields Tcs ~ 3 x 10^ 
K and L « 2600 Lq (iKahabka et al.|[l999i) . which IS con- 
sistent with a 0.5 Mq He WD near its maximum Toff. 
Neither of these sources seems to fit the canonical ex- 
planation of supersoft x -ray sources as th ermal timescale 
mass transfer systems ([King et al.ll200l"[) . Furthermore, 
IPietsch et all ((20051 l2007D 's growing catalog of M31 su- 
persoft sources may also prove to be useful in this search 
for He WD post-nova supersoft x-ray sources. 

6. CONCLUSIONS 

Prompted by the lack of convincing detections of He 
WDs in CVs, we have studied their fate and thermal 
evolution in the hopes that observers may find new tech- 
niques with which to reveal or constrain this population. 



We noted the existence of pre-CVs with He WDs in fj2] 
and derived the mass-transfer evolution for two examples 
in ^2.2l under the assumption of angular momentum loss 
due to gravitational wave radiation. In ^ we described 
the thermal evolution of the core both analytically and 
numerically for several cases with constant M, and in 
fJH we described the reheating histories for the time- 
dependent CV evolutions of ij2.2l In 53 we focused on 
the details of the individual nova events, describing the 
convective burning phase f i|5.1|) . the ejecta composition 
(t j5.2p . and the post-CN outburst quasi-steady burning 
phase f ^O)) . 

In fJH we addressed the issue of the discovery of these 
systems while in quiescence. We now turn to the ques- 
tion of their detection via their nova outbursts and post- 
nova burning phases . From populat i on synthesis predic - 
tions (Ide Koo]|ll99l IPolitanol[l996l ; IHowell et al.ll2001h . 
« 20% of CVs below the period gap harbor a He WD 
accretor. If, for simplicity, we consider only CNe on 0.4 
and 0.6 Mq WDs, whose ign ition masses differ by a fac- 
tor of « 2 f or a g iven Porb (jTownslev fc BildstenI 120051 ; 
lYaron et al.l[200l . we find that the rate of CNe from 
He WDs is « 10% of all CNe with orbital periods be- 
low the period gap. However, for some novae, and es- 
pecially extragalactic novae, the orbital period is much 
harder to measure, and so we cannot always divide the 
nova population into systems above and below th e pe- 
riod gap. From theory and observations ( Warner 2002| ; 
iTownslev fc BildstenI [2005t ). CNe below the period gap 
make up 10 — 25% of all CNe, and thus we expect He 
WD accretors in 1 — 3% of all CNe with no Porb informa- 
tion. Note that this rough estimate does not take into 
account the selection effects that arise from the different 
peak luminosities and outburst timescales, which depend 
on the CV parameters. 

As we showed in ^5.31 the He WD remains bright in UV 
and supersoft x-rays during the post-nova burning phase. 
Due to the similar physics at play (i.e., H-burning on 
low-mass degenerate cores), these post- novae have simi- 
lar observed characteristics to low-mass post-asymptotic 
giant branch ( AGB) and post-early AGB stars. From 
IBuzzoni et al.l (j2006l ) and references therein, the post- 
AGB birth rate is 2 per yr in a 10" Mq E/SO galaxy 
Extragalactic studies of CNe find a total rate of 20 ± 10 
per y r in a galaxy of the same size (IWi lliams fc Shafted 
I2004f ). of which, from the calculation above, 0.2 — 0.6 of 
these will be novae from He WDs. Thus, there should 
be 10 — 30% as many post-novae He WDs as post-AGB 
stars. 

As these post-novae cool, their evolution takes them 
towards the pulsational instability strip. lArras et al.l 
(|2006l ) examined the effect of envelope composition on 



He WDs IN CVs 



11 



pulsating WDs, and found that above a helium mass frac- We thank Dean Townsley for a thorough review of a 

tion of w 0.3, a Hell convection zone is present that re- preliminary draft, Dina Prialnik and Attay Kovetz for 
suits in a second instability strip, with a hotter blue edge use of their code, Boris Gansicke, Tom Marsh, and Ma- 
than the canonical H/Hel strip. If the post-nova envelope rina Orio for helpful discussions, and the referee for con- 
has a helium mass fraction of « 0.4 due to core-envelope structive comments. This work was supported by the 
mixing. Arras et al, ( 2006,) predict the blue edge of the National Science Foundation under grants PHY 05-51164 
resulting instability strip to be at « 1.8 x 10"* K, far and AST 07-07633. 
hotter than the standard empirically-derived strip. 

APPENDIX 

Given the long accumulation timescales between novae on low-mass He WDs, we provide here an analytic estimate 
of the importance of chemical diffusion from the solar composition envelope into the He core. If the diffusion timescale, 
Td, is of order or less than the accumulation timescale, chemical diffusion will be non-negligible. In the limit that 
is much shorter than the accumulation timescale, the composition will be in a diffusive equilibrium governed by a 

balance of the gravitational and electric forces. For a degenerate He gas, the electric field is eE = —2mpg, where g is 
the gravitational acceleration, and so the net force on a hydrogen nucleus is —iripg. In diffusive equilibrium, balancing 
the forces yields an exponentially decreasing hydrogen tail going into the He core with a length scale roughly equal to 
the pressure scale height. Heavier elements like He, ^^C, ^''N, and ^^O, on the other hand, experience no net force in 
the degenerate core because they have the same charge-to-mass ratio as helium, and so they perform a random walk 
inwards and will have a constant mass fraction in the core when in diffusive equilibrium. 

The timescale governing diffusion ov er a scale height is T d H '^/D, where H = P/pg is the pressure scale height, 
and D is the diffusion coefficient. From lAlcock fc lUarionovl (|1980l ). the diffusion coefficient for ions in a perfect gas of 
background plasma is 

3{2kTf/^ _ cm^ T^/^aY^ 

" 16ni(7rmi)i/2z2z2e4A, s p^Z^ZfA^ 

where = p/10'^ g cm^'^, T7 is the temperature in units of 10^ K, Z and A are the ion's charge and atomic weight, 
the Coulomb logarithm is A, the 1-subscript refers to the background ion, and the i-subscript refers to the diffusing 

ion species. For the case of a degenerate He core, the scale height is H — 6.2 x 10^ cm p^^^i^/gr), where 57 is the 
gravitational acceleration at the core-envelope interface in units of 10^ cm s^^. For a 10"'^ Mq envelope on a 0.4 Mq 
core, gr ~ 5. 

The effects of electron degeneracy are accounted for in the Coulomb logarithm. For electron degenerate material. 



A,; = In 



^ ^ (fcr)3,^V^exp(r;) 1 ^ Tfrj'^'eMv) 



PsZf 



(2) 



where ry is the degeneracy parameter and the second line has been evaluated for a background of helium. The Coulomb 
ratio of Coulomb to thermal energies, P^ = iJcouiombZ-Ethcrmai, for He-He coupling near the outer edge of the core is 

1 /3 

0.57/93 /'^7j SO the helium is roughly an ideal gas. With these approximations, the diffusion timescale through a scale 
height becomes 

...2.4xl0»,r,;«(l)'||. (3) 

For a 10"'^ Mq envelope halfway through its accumulation phase, the base temperature is 5 x 10^ K, the density at 
the base is 3 x 10^ g cm~^, and so the diffusion timescale for H into the He core is w 3 x 10^ yr, assuming 77 = 2. 
For carbon diffusing inwards, Td ~ 10^ yr. Thus, for accumulation timescales < 3 x 10^ yr, the composition profiles 
will not be in diffusive equilibrium. 

Figure [13] shows the hydrogen (top panel) and helium (bottom panel) mass fractions as a function of pressure at 
different times during the 10* yr accumulation phase of a 0.4 Mq He WD accreting at a rate of M — 10~^^ Mq yr~^. 
From lower to higher pressure, these profiles are at times of 10 , 3 x 10^, 5 x 10'', 7 x 10'', and 9 x 10'' yr after the 
previous CN. Because the accumulation timescale is 10* yr, neither the H nor heavier elements should be in diffusive 
equilibrium. Indeed, as Figure [13] shows, the H is not yet in diffusive equilibrium; its exponentially-decreasing tail 
extends only w 50% of a pressure scale height inwards. 
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